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ABSTRACT 

Context. The quasi-Hilda comets (QHCs), being in unstable 3:2 Jovian mean motion resonance, 
are considered a major cause of temporary satellite capture (TSC) by Jupiter. Though the QHCs 
may be escaped Hilda asteroids, their origin and nature have not yet been studied in sufficient 
detail. Of particular interest are long TSCs/orbiters. Orbiters - in which at least one full revolution 
about the planet is completed - are rare astronomical events; only four have been known to 
occur in the last several decades. Every case has been associated with a QHC: 82P/Gehrels 3; 
lllP/Helin-Roman-Crockett; P/1996 R2 (Lagerkvist); and the possibly QHC-derived D/1993 F2 
(Shoemaker-Levy 9, SL9). 

Aims. We focus on long TSC/orbiter events involving QHCs and Jupiter. Thus we survey the 

known QHCs, searching for further long TSCs/orbiters over the past century. 

Methods. First, we confirmed the long TSC/orbiter events of 82P, 11 IP, and 1996 R2 in order 

to test our method against previous work, applying a general ,/V-body Newtonian code. We then 

used the same procedure to survey the remaining known QHCs and search for long TSC/orbiter 

events. 

Results. We newly identified another long TSC/orbiter: 147P/Kushida-Muramatsu from 1949 
May 1 4^ jo jj^ s - 1 96 1 July 1 5 . Our result is verified by integrations of 243 cloned orbits which take 
account of the present orbital uncertainty of this comet. The event involves an Li — > L\ transition 
as with 82P and 1996 R2; this may represent a distinct subtype of TSCs from QHC derived 
(Lj — ») longer captures exemplified by 1 1 IP and (probably) SL9, though this classification is still 
only based on a small database of TSCs. 

Conclusions. This is the third long TSC and the fifth orbiter to be found, thus long TSC/orbiter 
events involving Jupiter have occurred once per decade. Two full revolutions about Jupiter were 
completed and the capture duration was 12.17+jj- 2 .' years; both these numbers rank 147P as third 
among long TSC/orbiter events, behind SL9 and 11 IP. This study also confirms the importance 
of the QHC region as a dynamical route into and out of Jovian TSC, via the Hill's sphere. 

Key words, minor planets, asteroids formation - comets: general - celestial mechanics 
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1. Introduction 



Among all the asteroids that have been recorded up to the present in the asteroid database, e.g., the 
"JPL Small-Body Database" (http://ssd.jpl-nasa.gov/), a large number (more than 1000 
including unnumbered objects) are known to populate the region of the 3:2 mean motion resonance 
(MMR) with Jupiter, in the outer main belt. These are the "Hilda asteroids" (Schubart [19681 119821 
[19911 Ip [T9761 Yoshikawa [19891 Franklin et al. [19931 Nesvorny & Ferraz-Mello [1997l l. Their 
semimajor axes, a, concentrate in the range 3.7 AU < a < 4.2 AU, at eccentricities e < 0.3, and 
inclinations i < 20° (Zellner et al. 119851 ). This results in a range of the Tisserand parameter with 
respect to Jupiter, 7), of ~ 2.90-3.05, where 7) = cij/a + 2 ^a/aj(l - e 2 )cos/, with cij being the 
semimajor axis of Jupiter and / the mutual inclination between the orbits. The critical arguments for 
the Hildas, <p = 3/lj - 2/1 - m, librate about 0°, being stable in the long-term (A is mean longitude, vj 
is longitude of perihelion, and J indicates Jupiter). As regards physical properties, the low-albedo 
D- and P-types are more abundant in Hildas' surface colours than the small fraction of C-types 
(Dahlgren & Lagerkvist [19951 Dahlgren et al. [T997l Gil-Hutton & Brunini I20081 Licandro et al. 
2008 ). The surface colour of D- and P-type asteroids, such as Hildas and Trojans in the outer main 
belt, corresponds well with that of cometary nuclei (Fitzsimmons et al. 1994; Jewitt l2002l . which 
means that they are covered with a similar mineralogical surface, suggestive of a common origin. 

More than 50 Jovian irregular satellites are known at present (Jewitt & Haghighipour 2007). 
Cuk & Burns (2004 ) pointed out that the progenitor of the main prograde cluster, the Himalia fam- 
ily, was plausibly derived from the Hildas long ago, if it was captured by a gas-drag assisted mech- 
anism. Thus tracing the origin and nature of irregular satellites is very significant for studying the 
accretion processes in the early solar system. Satellite capture mechanics in the circular restricted 
three-body problem (CR3BP) or the A^-body problem has often been investigated (e.g., Henon 
[19701 Huang & Innanen [T9831 Tanikawa [T9831 Murison [T9891 Brunini et al. |T996l Nesvorny et 
al. 2003 2007). Reflectance spectra of Jovian irregulars, being dominated by D- and C-types (Luu 
1991; Grav et al. 12003 ), are comparable to those of Hildas. 

During the past half century, several Jupiter family comets (JFCs; cf. Levison 1996 ) have stayed 
in or near the Hilda zone, although being in unstable 3:2 MMR with Jupiter. Some of them have 
been transferred from outside to inside, vice versa, or from inside to inside of Jupiter's orbit by 
undergoing a temporary satellite capture (TSC) by Jupiter (e.g., Carusi & Valsecchi |1979l Tancredi 
et al. [19901 . Such a JFC is called a "quasi-Hilda comet" (QHC) by Kresak (fT979l , who identified 
three such objects: 39P/Oterma, 74P/Smirnova-Chernykh, and 82P/Gehrels 3. Di Sisto et al. ( 120051 ) 
integrated the motions of 500 fictitious Hildas for ~ 10 9 years, and found that most of them escaped 
from the Hilda zone into the JFC population, i.e., left the 3:2 MMR and evolved quickly on to un- 
stable orbits: such a chaotic diffusion from the Hilda zone has also been demonstrated by Nesvorny 
& Ferraz-Mello ( 119971 ). In addition, large-scale collisional processes, such as the late heavy bom- 
bardment, might also release small bodies from the Hilda zone into JFCs, e.g., see Gil-Hutton & 
Brunini (2000). Hence, some QHCs may indeed be such escaped Hildas themselves. Recently, Toth 
(2006) updated the QHC list, finding a total of 17 members (see Section 3). This includes bodies 
such as the renowned Comet D/1993 F2 (Shoemaker-Levy 9, SL9) that have undergone TSC by 
Jupiter and then disappeared after colliding with the planet. Surface spectroscopic (or colorimetric) 
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measurements for QHCs have only been carried out for 82P (De Sanctis et al. 2000), the results of 
which also indicate a taxonomic D-type. 

Among the QHCs, 39P/Oterma was the first known to be temporarily captured by Jupiter, in 
1936-1938 (Marsden |19621 l. However, this comet flew through the region near Jupiter over a rather 
short time, during which the comet did not complete a full revolution orbiting about the planet. In 
contrast, unlike 39P's "fly-through" capture, there is a different kind of TSC, in which at least one 
full revolution about the planet is completed; we deal with these in the present paper. Following 
Kary & Dones (1 19961 1 we call such objects "orbiters". These are often characterized by a long 
capture with very small perijove distance, usually lasting for ~ 10 years or more. Not all orbiters 
become such long TSCs (> 10 yr), although of course they last longer than the fly-through type. 
SL9 is a representative case for both long TSCs and orbiters. This comet was pointed out to have 
possibly been QHC-derived before its tidal disruption on passing through perijove at less than 1 .5 
Jovian equatorial radii (Rj, where Rj - 71492.4 km), i.e., within the Roche limit, in 1992 July 
and its subsequent collision with Jupiter in 1994 July (Nakano & Marsden 1993; Sitarski 1995; 
Benner & McKinnon 1 19951 1. If it is QHC-derived, then SL9 is the only QHC so far that has been 
orbiting the planet as a TSC at the time of discovery (Benner & McKinnon 1 1 9951 l . By numerically 
integrating SL9's pre-collision orbital motion, several studies showed that the TSC duration of 
SL9 lasted for 50 years or more, during which the comet completed more than 30 revolutions 
orbiting about Jupiter, making it nominally the longest known TSC (Carusi et al. 1994; Benner & 
McKinnon 1 1995 1 Chodas & Yeomans 1996). However, according to Benner & McKinnon d 19951 ), 
SL9 was the most chaotic known object in the solar system with an effective Lyapunov time of 
only ~ 10 years on its jovicentric orbit. Thus it is difficult to determine with certainty SL9's pre- 
capture orbit and its true TSC duration. Nevertheless, Benner & McKinnon's additional statistical 
analysis of distributions in a-e space and Tj values, based on back integrations of the various SL9 
fragments, revealed a possible QHC origin of SL9. The work of Kary & Dones ( 1996 ) supports this 
possibility: they traced the motions of numerous fictitious JFCs for ~ 10 5 years and found that half 
the SL9-like very long captures > 50 years were due to QHCs. They estimated that the frequency 
of such a very long TSC (= "long capture" as designated by them) is extremely rare, only 0.02% 
of all the TSC events in their simulations. Interestingly, impacts on Jupiter are more frequent than 
the very long TSCs by a factor of 8-9. They also evaluated that long TSCs (= "orbiters bound > 10 
yr" as designated by them) and orbiters are still rare events, at the respective levels of 0.8% (d very 
long TSCs) and 2% (d long TSCs) relative to all events, with about 98% being the short TSC type 
which contains the fly-through events. Carusi & Valsecchi (119791 ) had earlier confirmed the rarity 
of orbiter events, simulating motions of fictitious small bodies as well. 

In the jovicentric Keplerian system, a TSC (especially a long TSC/orbiter) occurs whenever 
a small body passes near one of the collinear libration points L\ or L2 in the CR3BP of the Sun- 
Jupiter-(third) body system with very low velocity, i.e., effectively becoming bound by Jupiter 
when it enters the Hill's region with near-zero velocity. After that, the bound small body revolves 
about Jupiter on an elliptical jovicentric orbit until it again passes through the region near either 
L\ or L2 and escapes from the Jovian system. Considering such a transition using newly developed 
dynamical systems techniques based on a Hamiltonian formulation in the CR3BP, Koon et al. 
(2001) and Howell et al. (120011 ) demonstrated that a TSC by Jupiter occurs when the small body 
passes through a region inside the invariant manifold structure related to periodic halo orbits around 
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L\ or L2 in the Hill's region. The TSC (or its duration) is usually defined by the jovicentric Kepler 
energy, Ej, being negative, Ej < 0, with the additional condition that the bound small body must be 
within the jovicentric sphere of gravitational influence: Kary & Dones ( 1996 ) set its boundary at 3 
Hill's sphere radii (=1.065 AU) of Jupiter. However, Howell et al. (120011 1 defined the TSC duration 
as the residence time in the Hill's region. The former generally lasts longer than the latter, and here 
we regard the former as the TSC duration. The dynamics involved in TSC is quite different from 
that of quasi-satellites in 1 : 1 libration with Jupiter far outside the Hill's region (Wiegert et al. 2000; 
Kinoshita & Nakai|2(j07j). 

Apart from SL9, only three orbiters have been known to occur. Every case has been associated 
with a QHC: 82P; 1 1 lP/Helin-Roman-Crockett; and, though it is not in Toth's d2006b QHC list, 
P/1996 R2 (Lagerkvist). These TSCs were found by Rickman (119791 , Tancredi et al. ( fT990l ), and 
Hahn & Lagerkvist ( 1999 ), respectively, and are discussed further in Section 2. The occurrence of 
these few events during the last several decades is consistent with the rarity of long TSC/orbiter 
events suggested by Kary & Dones ( 119961 ). Although 74P encountered Jupiter at distances of 0.24 
AU and 0.47 AU in 1955 October and 1963 September respectively, it was not bound to Jupiter 
(Rickman 1979; Carusi et al. |1985bb ; this comet is, however, expected to experience a TSC by 
Jupiter in this century (Carusi et al. 1985b and see also Section 4). There further exist some JFCs 
that encountered Jupiter more closely than several QHCs involved in TSC, e.g., 16P/Brooks 2; 
D/1770 LI (Lexell); 8 IP/Wild 2; but they were not captured by the planet owing to their high- 
velocity encounters (Carusi et al. 1985a; Emel'yanenko 2003). 

Therefore studying the origin and nature of QHCs is of great importance from various as- 
tronomical points of view mentioned above, especially as regards their origin and being possibly 
related to the Hildas and the Jovian irregular satellites. Such studies may provide unique knowledge 
and clues about formation processes in the early solar system. Here we focus on long TSC/orbiter 
events involving Jupiter and QHCs. First (Section 2), applying a general A^-body Newtonian code, 
we reconfirmed the long TSC/orbiter events of 82P, 1 1 IP, and 1996 R2. Then (Section 3) we used 
the same procedure to search Toth's (2006) QHCs list for other objects that have become long 
TSCs/orbiters in the past century. Eventually, we successfully found another long TSC/orbiter, oc- 
curring in the mid-20th century, 147P/Kushida-Muramatsu. 

2. Computation method and its application to the known TSCs 

The jovicentric trajectories for the long TSCs/orbiters, except for SL9, have been numerically sim- 
ulated by Carusi & Valsecchi ( [1979] ) and Carusi et al. J1985a|[T985bl for 82P, by Tancredi et al. 
( [1990] ), Belbruno & Marsden ( fT997b and Howell et al. d200Tb for 1 1 IP, and by Hahn & Lagerkvist 
( [19991 for 1996 R2. 

First of all, we attempted to reproduce their TSC events in order to compare our simulations 
with these previous studies, integrating back to the time of each TSC. We applied a general A^-body 
Newtonian code, then equations of motion for the ith body are: 



where G is the gravitational constant, Mq and m, are respectively the mass of the sun and the ith 
body (i = 1 , 2, ■ ■ • , AO, r, is the heliocentric position vector of the ith body, r, is the heliocentric 
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distance of the z'th body, and r,j is the distance between the ith and j'th bodies. Here, we regarded 
comets as massless bodies. 

The integrator that we applied here is the "SOLEX", Ver. 9.1 package, developed by Vitagliano 
( 11997 !) based on the Bulirsh-Stoer method. Coordinates and velocities of the planets, regarded as 
point masses, were based on the JPL Planetary and Lunar Ephemeris DE409. We confirmed that the 
results of our numerical integrations did not significantly change when we used other integration 
methods that we have often applied in our studies, e.g., the Adams method. Our integrator can 
accurately process very close encounters by means of a routine that makes automatic time step 
adjustments, and truncation and round-off errors are almost negligible for our investigation here. 
Therefore the SOLEX integrator is sufficiently reliable to deal with the problem of close encounters 
with Jupiter. 

As initial parameters, up to date osculating orbital elements were taken from the JPL Small- 
Body Database, mentioned above in Section 1, for 82P and 1996 R2 and from Nakano (2005J for 
11 IP, as listed in Table Q] 82P and 11 IP are numbered, multiple-apparition comets, covered by 
very long observational arcs, hence their orbital solutions are highly precise. Meanwhile, although 
1996 R2 is unnumbered and still a one-apparition comet, consequently with a shorter arc, the num- 
ber of astrometric positions is the most, 135, and the RMS residual is the least, 0".76, among the 
three comets, hence we judged that the TSC of 1996 R2 is worth simulating here. The nongravita- 
tional parameters for these QHCs' motions were not detectable from their astrometry hence were 
not included in our integrations - in any case their motions are chaotic during TSC. The Jovian 
oblateness terms and the perturbations by the Galilean satellites were also ignored since their ef- 
fects are negligible in our investigations (see Kary & Pones fl996i l; these effects should be included 
in analyses of very long TSCs of SL9-like objects having highly inclined Jupiter-grazing orbits. 

Our simulated TSC trajectories of 82P, 11 IP, and 1996 R2 are shown in Figure Q] In these 
plots we use a rotating frame in which the jovicentric rectangular coordinates are referred to the 
heliocentric orbital plane of Jupiter, with the sun located in the direction of the -jt-axis, and the +z- 
axis pointing to the north. In Figure [1] every TSC motion is unstable and chaotic, with substantial 
solar perturbations occurring near each comet's apojove. 

Among these TSCs, the trajectories of 82P and 1996 R2 are similar in profile to each other. 
Both these QHCs came tangentially from outside Jupiter's orbit with a low-inclination retrograde 
motion, i.e., their pre-capture heliocentric orbits imply that they came from the Centaur region, 
their perihelia being beyond the Jovian orbit. During the TSC, first they passed close to L2 with 
low relative velocities (between Jupiter and comets in the jovicentric frame) of respectively ~ 0.7 
km s" 1 and ~ 0.6 km s _1 , then revolved about Jupiter, completing one full revolution, and finally 
escaped from the Jovian system into the QHC region, passing near L\. In the Lq — > L\ transition, 
an object that enters the Jovian system close to Li seems to have enough orbital energy to rapidly 
escape through the region near L\ (Tanikawa |1983l ). The long lasting captures for 82P and 1996 
R2 may be somewhat related to the orbital stability for the retrograde satellites (Huang & Innanen 
119831 Nesvorny et al. 120031) . However, the TSC duration is below 10 years in both cases; so to be 
more exact these orbiters should perhaps not be categorized as long TSCs. 

On the other hand, the motion of 1 1 IP was more complicated than those of 82P and 1996 R2. 
11 IP completed three full revolutions about Jupiter during its long TSC of 18.45 years, both of 
which numbers rank second next to SL9 among known TSCs. Its pre-capture and post-escape he- 
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liocentric orbits were both QHC type (Belbruno & Marsden |1997l l, entering through the region near 
L[ at a relative velocity of ~ 1.0 km s _1 and escaping near L\ again. Its capture/escape trajectory 
is almost symmetric about the jc-axis (cf. Murison 1989) inside the Hill's region in the CR3BP. 
Howell et al. (2001) demonstrated that the 11 IP capture/escape trajectory is almost identical to 
some asymptotic trajectory winding onto a periodic halo orbit of L\, such trajectories forming 
segments of the invariant manifolds. 

We therefore reconfirmed the long TSC/orbiter events involving 82P, 11 IP and 1996 R2. Our 
orbiter trajectories look exactly as in Belbruno & Marsden (1997) and Howell et al. (2001) for 
1 1 IP, and as in all the above papers for 82P, but the results by Tancredi et al. ( 1990) for 1 1 IP and 
Hahn & Lagerkvist (1999) for 1996 R2 show somewhat different behaviours. Both these can be 
ascribed to slight differences in initial parameters - their initial orbits were based on observational 
arcs of respectively 2—4 months and 111 days whereas ours were based on longer time spans of 
4643 and 153 days. Besides, Hahn & Lagerkvist (119991 ) plotted the trajectory not of the nominal 
orbit but of a clone, i.e., a possible TSC trajectory within the uncertainty of the initial parameters. 
In fact, we could not reach the TSC of 1996 R2 by using their nominal orbit. This demonstrates 
that in studying TSC events, initial parameters must be determined precisely, based on long-arc 
improved orbital solutions wherever possible, since jovicentric orbits during TSC have a chaotic 
nature. 

The capture situations of these QHCs are summarized in Table [2] During these TSCs, 82P, 
1 1 IP and 1996 R2 experienced 2, 6 and 3 perijove passages, respectively. In particular, 82P passed 
its closest perijove of 3.01/?j in 1970 August, just outside the Roche limit for comets ~ 2.1R] 
(assuming a cometary bulk density p ~ 1 g cm -3 ). Table|2]also shows the minimum values of the 
jovicentric Kepler energy Ej - -l/2a' (in ALT 1 ) where a' is the semimajor axis for the jovicentric 
orbit. The low values attained, -2.80, -3.20, and -3.28 for the three comets, are often a feature of 
long TSCs/orbiters, as against -1 .79 which we have computed for 39P during its fly-through in the 
1930s. 

Shortly after TSC, these three QHCs were discovered and observed as "quasi-Hilda comets" for 
the first time. The time difference, A^sc-obs)) between the end time of the TSC, given in Table [2] 
and the first astrometric time, given in Table [T] is notable. Every AT^sc-obs) is short enough (1.3 
years for 82P, 3.49 years for 1 1 IP, and 3.16 years for 1996 R2) to indicate that the simulated TSC 
should be a real event. 

3. 147P/Kushida-Muramatsu: another long TSC/orbiter 

Next we surveyed the remaining objects in Toth's (2006) QHC list (Tabled to determine whether 
other long TSCs/orbiters around Jupiter exist. We applied our code to integrate each orbit back 
100 years from its initial epoch, since the Lyapunov time of JFCs is rather short, generally 100 + 
50 yr (Tancredi I1995I ). Moreover, going back more than 100 yr may not be meaningful because 
of likely weak nongravitational accelerations, in addition to the accumulation of orbital errors. We 
recognized that short TSC events have sometimes occurred, for which results will be published 
as an additional paper elsewhere. Eventually, we successfully found another long TSC, involving 
1 47P/Kushida-Muramatsu . 
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The observational history of 147P is described by Green JT9931 |2000j i and Kronk ( |2008l ). The 
comet was photographically discovered by Y. Kushida and O. Muramatsu ( Yatsugatake South Base 
Observatory, Japan) on 1993 December 8.65, using a 0.25-m f/3.4 reflector. It was then at magni- 
tude ~ 16.5 and slightly diffuse, about 1-2' in diameter with a central condensation. The comet was 
recovered on 2000 October 3.72, when T. Oribe (Saji Observatory, Japan) obtained CCD images 
with the 1.03-m reflector. Thus 147P has been detected at two perihelion returns, allowing it to be 
numbered, although it has not been astrometrically observed since 2002 March. 

In our integration of 147P, the initial parameters were taken from Nakano (120021 1. as listed in 
Table |4] His orbital solution was accurately improved based on 197 carefully selected positions 
covering an arc of 3024 days and also provided associated error estimates (Table HJ. We believe 
that this solution is the best determined for this comet, probably based on a longer arc and more 
astrometric positions than any other. The nongravitational accelerations in 147P's motion were not 
detected from the astrometry. Integrating back this orbit, we found that the comet was captured by 
Jupiter in the mid-20th century. The 3D-view of 147P's TSC situation is illustrated in Figure|2] 

However, AT^sc-obs) of 147P is extremely large, more than 30 years, which greatly exceeds 
the values of < 3.5 years for 82P, 11 IP and 1996 R2 (and also SL9 with negative A!T(Tsc-obs))- 
Thus we must further examine whether or not 147P really underwent a long TSC/orbiter phase in 
the mid-20th century. To answer this, we traced other possible orbital motions of 147P back to and 
during its TSC, generating multiple "clones" at initial epoch, and integrating them. As well as the 
nominal 147P osculating orbit, other clones had slightly different orbital elements, within the +lcr 
error, given in Table |4] We generated the clones on the basis of three possible values [nominal and 
+ lcr] for five orbital elements [a; e; a>; Q\ z']. This number of permutations yielded a total of 243 
(= 3 5 ) clones, including the nominal one. We also generated some test clones varying T (or M) 
within the allowed error, confirming that this effect is negligible; thus we excluded this parameter 
when generating the main set of clones. Conversely, we found a to be the most strongly effective 
parameter in causing divergence among the clones' orbital behaviours. 

We found that all 243 clones underwent TSCs. In Figure [3] we can see that their trajectories 
become scattered in their pre-capture orbital phase, while converging toward the post-escape phase, 
since the orbital motions were integrated back in time. However, the dispersion of their TSC tra- 
jectories, corresponding to +lcr error in the orbit determination, is not as scattered as is sometimes 
expected with chaotic trajectories. We conclude that it is very likely that 147P really was an or- 
biter, experiencing a long TSC in the mid-20th century. The motions of all the clones follow almost 
the same profile as those of 82P and 1996 R2, coming tangentially from outside Jupiter's orbit, 
i.e., with a low-inclination retrograde motion in the jovicentric frame. The pre-capture heliocentric 
Centaur orbit of 147P had a ~ 6.2 AU, q ~ 5.8 AU, and e ~ 0.07. The comet entered the Jovian 
system passing through the region near Li with relative velocity ~ 0.9 km s _1 , revolved about 
Jupiter, escaped from the Jovian system passing near L\, and was discovered and observed as a 
QHC after ~ 32 years. However, an important difference from both 82P and 1996 R2 is that 147P 
completed two full revolutions about Jupiter, thus experiencing a significantly longer TSC. 

Figure |4] shows the variations in rj,£j, Eq, and Tj of all the clones in and around the TSC 
interval, where rj is jovicentric distance (in AU), and Eq is heliocentric Kepler energy = -1 /2a (in 
AU -1 ). The rj and Ej diagrams in the pre-capture orbital phase display a scatter that is as evident as 
the scatter in Figure [3] We can see (most clearly from the rj diagram, but with signs also reflected 
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in the other three plots) that 147P experienced three perijove passages. The term of negative Ej 
corresponds to TSC, during which rj was always within 3 Hill's sphere radii, as defined in Section 
1; thus 147P underwent a long TSC for H.ll^^ years (where error estimates are based on the 
clones' dispersion). The minimum Ej taking a low value of -2.78 + 0.01 is also an expected feature 
of a long TSC/orbiter event. In the Eq diagram, we can see the sharp spike corresponding to the 
closest encounter to Jupiter, during which the nominal heliocentric orbit is briefly hyperbolic. The 
value of 7j was slightly higher than 3.0 in the pre-capture orbital phase, as is often associated 
with low-velocity encounters at Jupiter (Kresak 1979), and it suddenly dropped twice around the 
closest perijove passages. We summarize the capture situation data of 147P in Table[5](again, error 
estimates come from the range shown by all the clones). 

Also indicated in the Eq diagram in Figure [4] is the commensurability of 147P with Jupiter 
changing from 3:4 Centaur-type in the pre-capture orbital phase to 8:5 (rather than 3:2) QHC in 
the post-escape phase. This, however, corresponds to the osculating orbit immediately post-escape, 
and integrations for motions of all the clones over a longer timescale suggest instead that 147P 
might librate in the 3:2 Jovian MMR for ~ 350 yr following the TSC (and is doing so at the present 
time), though with a somewhat faster libration period and higher libration amplitude (~ 135° for 
the nominal orbit) than the more stable, typical Hilda asteroid orbits (e.g., Schubart 1982, 1199T1 
Franklin et al. 1993; Nesvorny & Ferraz-Mello 11 9971 ) . Moreover, test integrations of the other 
objects suggest that of the QHCs known to have been involved in a TSC, only 147P has undergone 
an interval where the critical argument of the 3:2 MMR librates. 

Further backward integrations of 147P inform us of the possibility of another TSC somewhere 
between the late- 19th to early-20th century. However, as all the clones were extremely scattered by 
that time, we do not consider that TSC here. 

4. Concluding Remarks 

On the basis of our investigations above, we have presented a newly identified TSC of a comet by 
Jupiter. This is in the rare, orbiter class of TSCs and involves 147P from 1949 May 14^^^-1961 
July 15. This is the third long TSC of > 10 years and the fifth orbiter found, so that Jupiter's long 
TSCs/orbiters have occurred once per decade. The completion of two full revolutions about Jupiter 
and the capture duration of 12.17*0^ years rank 147P as third in both these numbers among known 
orbiters, behind SL9 and 1 1 IP. 

4.1. TSC Classification 

Following Kary & Dones (I1996l l, we classify the known TSCs as follows: 39P as fly-through, 82P 
and 1996 R2 as orbiters, 1 1 IP and 147P as long TSCs (> 10 years), and SL9 as a very long TSC 
(> 50 years). 

On the other hand, depending on the TSC characteristics, Howell et al. (120011 1 defined two 
TSC types: Type 1 as a 39P-like fly-through; Type 2 as a long lasting capture like 1 1 IP, where the 
comet experiences more than one close encounter with Jupiter while in the TSC region. Here we 
tentatively divide Type 2 into two subtypes: Type 2A as the L2 — > L\ transition such as 82P, 1996 
R2 and 147P; Type 2B as the QHC-derived (Li -->) longer capture than 2A, as with 1 1 IP and SL9. 
Further TSC classifications may be possible in the future, if different kinds of TSC are found. 



Ohtsuka et al.: Quasi-Hilda comet 147P/Kushida-Muramatsu 

4.2. Future TSCs 
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We also surveyed all the known QHCs, integrating their orbital motions forward for 100 years to 
check for future long TSCs/orbiters. We found that 11 IP will undergo a long TSC/orbiter phase 
with 6 perijove passages, and minimum Ej ~ -3.19, from 2068 April 20-2086 June 09 (duration 
~ 18.14 years). We identified a rather long capture of 82P with minimum Ej ~ -2.26 though it 
is not an orbiter but instead follows a symmetric trajectory about the x-axis, from 2056 February 
18-2064 July 26 (duration ~ 8.43 years). There appear to be two TSCs for 74P if we take its initial 
parameters from the JPL Small-Body Database, from 2025 August 10-2031 May 29 (duration 
~ 5.80 years, minimum Ej ~ -0.74) and 2081 January 29-2085 January 29 (duration ~ 4.00 
years, minimum Ej ~ -0.77). 

Such future TSCs have already been predicted by Carusi et al. dl985a|[T985bt for 74P and 82P 
and Tancredi et al. ( 1990J for 1 1 IP, and our simulations are almost identical in profile with those. 
Meanwhile, it is unlikely that 1996 R2 and 147P will encounter Jupiter as a long TSC/orbiter for 
the next 100 years from each initial epoch. No further long TSCs/orbiters were found among the 
remaining QHCs. 

4.3. Tidal Effects 

Here we discuss the Jovian tidal force acting on 147P and other captured QHCs, around their 
perijove passages. Well-known cometary tidal splitting events have occurred twice: 16P in 1886 
(Sekanina & Yeomans 1985) and SL9, both passing perijove within the Roche limit for comets, 
~ 2.7 Rj. In contrast, the close encounter of 147P with Jupiter around 1952 August 26 (Table [5]) 
was at a perijove distance of 14.61*?'?ji?j, comfortably outside the Roche limit. The radius, R c , 
of 147P's nucleus was estimated at 0.21 km by Lamy et al. (f2()04), the smallest among all their 
measured cometary nuclei. 

The Jovian tidal stress, cr T , acting on small bodies (QHCs here), is given by: 

ctt * Gm sP R 2 Jrl (2) 

where ;«j is the mass of Jupiter. Compared to the <x T acting on SL9 at perijove < l.5Rj in 1992 
(Scotti & Melosh [19931 Asphaug & Benz [19941 [19961 ), the cr T for 147P was extremely small, 
amounting to < 0.005% if we take the original, pre-encounter R c of SL9 as ~ 1 km (Scotti & 
Melosh [19931 Asphaug & Benz 1994) and p of both the comets as being equal. Hence, although 
our treatment here is approximate, there seems no reason to believe that 147P was affected by 
Jovian tides. On the other hand, 82P, which passed its perijove at 3.01/?j in 1970, just outside the 
cometary Roche limit, may have suffered some effect due to the Jovian tides, even though it was 
not broken up. Taking R c for 82P as 0.73 km (Lamy et al. 120041 ). and equal p as above, implies 
that the cr T acting on 82P could be up to 7% of that for SL9. If 82P is a friable and furthermore a 
loose rubble-pile object, the heating energy from continual Jovian tides might sufficiently affect the 
comet's nucleus structure so as to allow H2O ice, if present, to sublimate even though 82P is beyond 
the usual heliocentric distance for H2O sublimation. If our hypothetical scenario is true, such tidal 
heating effects could trigger cometary activity on 82P and produce an outburst. In this event, the 
coma size and consequent brightness would peak several days (or more) after the maximum tidal 
effects. 
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Another intriguing point is whether or not other tidal effects, e.g., tidal distortion and tidal 
torques leading to rotation state changes (Scheeres et al. 2004), are detectable in the light-curve 
observations of the captured QHCs. In the physical database of cometary nuclei by Lamy et al. 
(2004), we can notice a rather high axis ratio (a/b > 1.6) and long rotation period ~ 50 hr in the 
dataset of 82P, though an interpretation in terms of such tidal effects is still speculative. 

4.4. Comets or asteroids? 

At any rate, we have considered here only a small subset of the TSC events, focusing on the long 
TSCs/orbiters, for which there still exist only limited available data. Therefore we know little about 
such unique and extraordinary astronomical events, which are still full of ambiguities. In addition, 
we do not know whether QHCs are comets or asteroids. Indeed, several QHCs have sometimes been 
discovered or recovered as asteroids because of their cometary activity being weak: e.g., 36P (= 
1925 QD = 1940 RP), 39P (= 1950 CR), 74P (= 1967 EU = 1978 NA 6 = 1981 UHi 8 = 1982 YG 3 ), 
D/1977 CI (= 1977 DV 3 ), P/1999 XN 120 , P/2001 YXi 27 , and P/2003 CP 7 . Further observations 
and research for the QHCs will be necessary to unlock their origin and nature. 
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Table 1. Initial parameters of QHCs 82P/Gehrels 3, 1 1 lP/Helin-Roman-Crockett, and P/1996 R2 
Lagerkvist (equinox J2000) 



object 



82P 



111P 



1996 R2 



osculation epoch (TT) 
mean anomaly M 
semimajor axis a (AU) 
eccentricity e 
argument of perihelion w 
longitude of ascending node Q 
inclination i 

number of astrometric positions 
astrometric arc 

RMS residual 
source reference 



2002 Mar 27.0 
24°. 27332 
4.1388678 
0.1238743 
227°.65115 
239°. 62851 
1°. 12654 
103 

1975 Oct 27-2002 Mar 08 
9629 days (26.36 yr) 
0".79 
JPL 12 



2004 Dec 21.0 
359°. 25243 
4.0400597 
0.1402634 
10°.56557 
91°.93769 
4°. 23300 
126 

1989 Jan 03-2001 Sep 20 
4643 days (12.71 yr) 
0".86 
Nakano (2005) 



1996 Oct 08.0 
346°. 17074 
3.7820178 
0.3099920 
334°. 045 13 
40°. 24084 
2°. 60532 
135 

1996 Aug 12-1997 Jan 12 
153 days (0.42 yr) 
0".76 
JPL 18 
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Table 2. Long TSC/orbiter data of QHCs 82P, 1 1 IP, and 1996 R2 



object 


TSC duration 


no. 


min. 


heliocentric orbit & 


perijove time 


dist. e 


i / f 
Vj 




(TT) 


rev." 


E, b 


Li to Lj transition 


(TT) 


(AU) 


(km s" 1 ) 


82P 


1966 Dec 11-1974 Jul 11 


1 


-2.80 


2 : 3 -> 3 : 2 


1970 Aug 15.6 


0.0014 (3.01 Rj) 


34.25 




7.58 yr 








1973 Mar 24.9 


0.041 (85.16 if,) 


9.18 


1 1 IP 


1967 Jan 25-1985 Jul 08 


3 


-3.20 


3:2^3:2 


1969 Oct 14.4 


0.057 (119.4 Rj) 


5.00 




18.45 yr 






U ->Li 


1972 Sep 03.2 


0.175 (336.4 Rj) 


2.33 












1976 Apr 11.2 


0.012(24.31 R,) 


11.85 












1979 Jul 18.6 


0.211 (442.4 Rj) 


2.11 












1980 Nov 17.0 


0.201 (421.1 R,) 


2.20 












1983 Aug 10.1 


0.063 (131.7/?,) 


4.74 


1996 R2 


1983 Sep 07-1993 Jun 13 


1 


-3.28 


7 : 13 -> 8 : 5 


1987 Mar 19.7 


0.0074 (15.53 Rj) 


14.92 




9.77 yr 






L 2 -» L, 


1990 Jun 28.8 


0.169 (352.7 Rj) 


2.78 












1992 Sep 18.5 


0.325 (680.0 R s ) 


1.82 



a number of completed full revolutions about Jupiter 

b minimum E s in AIT 1 - occurred on 1970 July 30 for 82P, 1976 February 11 for 11 IP, and 1987 March 
15 for 1996 R2 

c Transition associated with the TSC: firstly, the Jovian MMR (heliocentric orbit) the comet is closest to 
immediately before/after the TSC; and secondly, the L, point near which the comet passes at the start/end 
of the TSC phase. 

d time of perijove passage in TT 

c perijove distance in AU and also in Rj 

f jovicentric velocity in km s -1 at perijove 
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Table 3. Quasi-Hilda comets 



Numbered QHCs: 
36P/Whipple 
82P/Gehrels 3* 
135P/Shoemaker-Levy 8 
Unnumbered QHCs: 
D/1977 CI (Skiff-Kosai) 
P/2001 YX 127 (LINEAR) 



39P/Oterma* 

1 1 lP/Helin-Roman-Crockett* 
1 47P/Kushida-Muramatsu 

D/1993 F2 (Shoemaker-Levy 9)* 
P/2003 CP 7 (LINEAR-NEAT) 



74P/Smirnova-Chernykh 
117P/Helin-Roman-Alu 1 



P/1996 R2 (Lagerkvistf* 
P/2002 08 (NEAT) 



77P/Longmore 
129P/Shoemaker-Levy 3 



P/1999XN 120 (Catalina) 
P/2004 F3 (NEAT) 



a not included in Toth's f2006l > QHC list 
* involved in previously known TSC 
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Table 4. Initial parameters of QHC 147P/Kushida-Muramatsu (equinox J2000) and their + lcr error 
estimates (Nakano 2002) 



osculation epoch (TT) 
perihelion time T (TT) 

perihelion distance q (AU) 

semimajor axis a (AU) 

eccentricity e 

argument of perihelion a> 

longitude of ascending node Q 

inclination 

number of astrometric positions 
astrometric arc 

RMS residual 



2001 May 11.0 

2001 Apr 29.48509 ±0.00081 
(=> M = 1°. 52643) 

2.7524408 ±0.0000019 

3.8094146 ±0.0000010 

0.2774636 ±0.0000005 

347°.55482 ±0°.00047 

93°. 69336 ±0° .00044 

2°.36694 ±0°.00004 
197 

1993 Dec 08-2002 Mar 20 
3024 days (8.28 yr) 
0".85 
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Table 5. Long TSC data of QHC 147P 



TSC duration 
(TT) 


no. 
rev. 


min. 


heliocentric orbit & 
L t to Lj transition 


perijove time 
(TT) 


dist. 

(AU) 


Vj 

(kms- 1 ) 


1949 May 14-1961 Jul 15 
(= JDT 2433050.5^- 
JDT 2437495.5) 


2 


-2.78 
+0.01 


3:4^8:5 


1952 Aug 26.6^| 
1955 Jan 18.8^ 8 
1960 Jan 20.86±0.02 


0.0070~ (14.61+^j) 
0.0273^»JJ» (57.05!!^) 
0.220 (460.60±0.11/?j) 


15.45^ 

7 cq+0.22 
J -0.20 

2.73 



12.17^ yr 



occurred on 1955 May 28.0 (JDT 2435255.5) +6 days 
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I 1 1 I 1 1 I — I — 1 -0.5 

-0.5 0.5-0.5 0.5-0.5 0.5 

X (AU) 

Fig. 1. Simulated TSC trajectories for 82P, 1 1 IP, and 1996 R2. Plots show x-y projections, over 1 
AU square, onto the heliocentric orbital plane of Jupiter, in the jovicentric rotating frame where the 
sun is always in the direction of the -x-axis. The arrows indicate the direction of time (opposite to 
the direction of the integrations, which went backwards in time). 




Fig. 2. 3D-view of 147P's trajectory during TSC in the jovicentric rotating frame covering 1 AU 
square. The vertical lines, spaced at 10-day intervals, each connect a position of the comet on the 
trajectory with its projection onto the x-y plane. 
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Fig. 3. x-y and x-z projections of TSC trajectories for 243 clones generated from 147P. They un- 
derwent long TSCs, completing two full revolutions about Jupiter during the Li — > L\ transition. 
The dispersion of their trajectories corresponds to +l<x error in the orbit determination. Hence we 
conclude that the real 147P indeed experienced a long TSC/orbiter event in the mid-20th century. 
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Fig. 4. Variations of rj, £j, Eq and Tj of 243 clones of 147P in and around their interval of TSC, 
where rj = jovicentric distance (in AU); Ej and £0 = jovicentric and heliocentric Kepler energy, 
respectively; 7j = Tisserand parameter. The abscissa for all the plots covers 2431596.5 (1945 May 
21.0 TT) to 2438095.5 (1963 March 07.0 TT), the time for which the nominal 147P was within 2 
AU of Jupiter. 



